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Abstract 

Let D be a domain in the plane containing the origin. We are interested in the 
ensemble of self-avoiding walks (SAW's) in D which start at the origin and end on the 
boundary of the domain. We introduce an ensemble of SAWs that we expect to have the 
same scaling limit. The advantage of our ensemble is that it can be simulated using the 
pivot algorithm. Our ensemble makes it possible to accurately study SLE predictions for 
the SAW in bounded simply connected domains. One such prediction is the distribution 
along the boundary of the endpoint of the SAW. We use the pivot algorithm to simulate 
our ensemble and study this density. In particular the lattice effects in this density that 
persist in the scaling limit are seen to be given by a purely local function. 

1 Introduction 

In two dimensions there has been a lot of interest in the self-avoiding walk (SAW) in simply 
connected domains D because of its conjectured relationship with SLE. One is interested in 
two cases: the radial case in which the SAW starts at a point inside the domain and ends on 
the boundary and the chordal case in which it starts and ends at boundary points. If one fixes 
the starting and ending points then Lawler, Schramm and Werner conjectured that the scaling 
limit is radial or chordal SLE 8 / 3 [7]. If one fixes the starting point but allows all SAW's that 
end anywhere on the boundary, then there are conjectures for the hitting density along the 
boundary from SLE partition functions [TJ Ej. Some progress towards proving the conformal 
invariance of the SAW was made by Duminil-Copin and Smirnov [2]. 

For the SAW in the full plane or a half plane there is a very fast Monte Carlo algorithm 
known as the pivot algorithm. (Clisby's recent implementation of the algorithm using binary 
trees has dramatically increased its speed [I].) However the pivot algorithm cannot be used for 
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the SAW in a bounded domain. The pivoting will almost always produce a SAW that does not 
end on the boundary of the domain or leaves the domain. And the pivot algorithm works on 
an ensemble of SAWs with a fixed number of steps. For the SAW in D we must allow SAWs 
with all numbers of steps. 

In this paper we show how to use the pivot algorithm to simulate an ensemble of SAWs that 
should have the same scaling limit as the ensemble of SAWs in D that end on the boundary. 
The key idea is to consider the ensemble of all SAWs starting at the origin such that there 
is a positive constant A > such that the SAW stays in the dilated domain XD and ends on 
the boundary of XD. We refer to this ensemble as the "dilation ensemble." The SAWs in the 
dilation ensemble have an arbitrary number of steps. Nonetheless we will argue that in the 
scaling limit this ensemble can be realized using the ensemble of walks with a fixed number of 
steps that satisfy the constraint above. This allows us to use the pivot algorithm to simulate 
the scaling limit of the ensemble. 

If one allows the SAW to end anywhere on the boundary of the domain, then in the scaling 
limit the distribution of the endpoint along the boundary gives a probability measure on the 
boundary analogous to harmonic measure for the ordinary random walk. Assuming the bound- 
ary is smooth, the distribution should be absolutely continuous with respect to arc length along 
the boundary. Lawler, Schramm and Werner have given an explicit conjecture for this boundary 
density [7] using SLE partition functions. The simplest description of their conjecture is that 
the boundary density for the SAW is proportional to the density for harmonic measure raised 
to the 5/8 power. In [7] they considered domains in which the sides are parallel to a lattice 
direction. In the case of a horizontal strip with the SAW starting on one boundary of the strip 
and ending on the other boundary, the conjecture for the boundary density was tested in [3], 
and good agreement was found. The conjecture was discussed further by Lawler in [6] where it 
was stated for general simply connected domains "after taking care of the local lattice effects." 

In [5] it was conjectured that the local lattice effect at a point z on the boundary only 
depends on the angle of the tangent to the boundary at z. There are a variety of ways to 
interpret what it means for the SAW which is on a lattice to end on the boundary of D which 
typically does not pass through lattice sites. The precise nature of the lattice effects depends 
on which interpretation is used. Explicit conjectures for these local lattice effects were given 
in [5] for two particular interpretations of ending on the boundary. In one interpretation one 
considers all SAWs u with u(i) G D for i < \u\ and co(|a;|) ^ D, where \u\ is the number 
of steps in the SAW. While this is one of the most natural interpretations, there is no good 
algorithm for simulating this ensemble. The other ensemble considered in [5] is the "cut-curve" 
ensemble. It uses infinite length SAWs in the full plane and conditions on the event that the 
SAW crosses the boundary of D only once. (In practice we use SAWs with a fixed length N 
and take N large enough that the typical size of the SAW is much larger than the domain D.) 
SAWs in this ensemble are equivalent to the concatenation of a SAW in D from the origin to 
the boundary of D and a SAW in the exterior of D, starting at the same boundary point and 
going to oo. 

One can use the pivot algorithm to simulate the "cut-curve" ensemble. But this has the 
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disadvantage that one must deal with a double limit - one must first let the number of steps in 
the SAW go to infinity and then take the lattice spacing to zero. In practice this means that 
one must simulate SAW's that are much larger than the domain D. So such simulations are 
less efficient than the method we will introduce for the dilation ensemble. 

In the next section we define the dilation ensemble and show how it may be simulated using 
the pivot algorithm on the ensemble of walks with a fixed number of steps. In section three we 
give an explicit conjecture for the lattice effects in our dilation ensemble and show how these 
effects may be computed by simulation. In section four we use our method for simulating the 
dilation ensemble to study the SAW in several different domains. We compare the boundary 
densities we find with the SLE partition function predictions for these densities. 

2 The dilation ensemble 

We first consider the radial case. We assume that our simply connected domain D contains the 
origin and has a boundary which is a smooth curve which we denote by dD. We also assume 
that the domain is star-shaped with respect to the origin, i.e., any ray from the origin intersects 
dD in only one point. (It would be interesting to generalize our approach to domains that do 
not have this property.) We are interested in the ensemble of SAWs that start at the origin 
and end on dD. Note that we are not fixing the point on the boundary where the SAW must 
end; the endpoint of the SAW is a random point on the boundary. 

The dilation ensemble can be described as all SAWs such that the domain D can be dilated 
so that the SAW lies in the dilated domain and ends on the boundary of the dilated domain. 
We let 5 > denote the lattice spacing. We start with the ensemble of all SAWs in the full 
plane on the lattice <5Z 2 that start at the origin with any finite number of steps. For a SAW uo 
we let \ui\ denote the number of steps in the walk. The walk co is given the weight where 
fi is the connective constant. (It is given by fi = limTv-i-oo cjx where Cat denotes the number of 
SAWs in the full plane with N steps that start at 0. The existence of this limit follows from 
the subadditivity of ln(cjv) |8J.) We will impose two constraints on the ensemble of all finite 
length SAWs. We let X(u) > be the dilation such that the endpoint of co is on the dilated 
curve \(u))dD. The partition function for the ensemble of all finite length SAWs starting at the 
origin is infinite, so we add the constraint that X(lj) lies in [Ai, A 2 ] where Ai < A 2 are positive 
constants. We expect that this makes the partition function finite, although we cannot prove 
this. The walk u ends on the curve X(cj)dD, but it need not stay inside the region X(co)D. Our 
second constraint is that the walk uj does stay strictly inside this region except for the endpoint. 
We let 1d(o;) denote the indicator function that is 1 if cu stays strictly inside this region except 
for its endpoint and if it does not. So our ensemble corresponds to the partition function 

Z = ^-Ml^w) l(Ai < A(w) < A 2 ) W{u) (1) 

We have included a weighting factor W(u>) which we will define later to make our ensemble 
correspond to the ensemble of SAWs in D that end on the boundary. Associated with this 
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partition function is a probability measure on simple curves in D that end on the boundary. 
Let ^(7) be a function on simple curves 7 in D that start at the origin and end on dD. If 
1d{w) = 1, then -^y is such a curve and so V'(a^)) * s defined. Define 

Zty) = < A(w) < A 2 ) W(u) ^(^-) (2) 

The expected value of ip is defined to be Z{ip)/Z. (Both Z and Z{ifj) depend on Ai and A 2 , but 
they will be fixed throughout our discussion.) 

We now consider how this ensemble corresponds to the ensemble of SAW's in D that end 
on dD. For a large integer N, let dX = (A 2 — Xi)/N. We think of the region between the 
curves X\dD and X 2 dD as the union of the regions between the curves (Ai + (k — l)d\)dD and 
(Ai + kd\)dD where k = 1, 2, • • • , N. In the scaling limit, the ensembles corresponding to these 
different regions are related by just a dilation. So they will all give the same expected value 
to if). The region between (Ai + (k — l)dX)dD and (Ai + kdX)dD corresponds to thickening 
the boundary of (Ai + kdX)dD. So we can think of the dilation ensemble as interpreting the 
constraint that the SAW stays in D and ends on dD as thickening the boundary. 

It is important to observe that the way that our ensemble thickens the boundary is not 
uniform along dD. The natural way to thicken the curve would be to take the thickness in 
the direction perpendicular to the curve to be constant along the curve. Let p(z) denote the 
density of the endpoint of the SAW in D which ends on the boundary. In the dilation ensemble 
if we take the weight W(u) to just be constant, then the corresponding boundary density will 
be proportional to p(z) times the thickness at z. To correct for this we take W(u) to be 
proportional to the inverse of this thickness. 

Recall that we assume that our domain is star-shaped with respect to the origin, i.e., any ray 
from the origin intersects dD in exactly one point. So we can parametrize dD by the polar angle 
9 with respect to the origin. We let D(9) be the distance from the origin to the curve in the 
direction 9. Along the ray at polar angle 9 the thickness is (A + dX)D{9) — XD{9) = D{9)dX. 
But this segment is not perpendicular to the tangent line. Let a(9) be the angle of a line 
perpendicular to dD at our point. Then the thickness of our shell perpendicular to the tangent 
line is D{9) cos ((9 — a{6)). Thus we define 

W{u) = [D{9{u)) cos(fl(w) - a(9(u)))]- 1 (3) 

where 9(u) is the polar angle of the endpoint of u. Note that the weight only depends on the 
polar angle of the endpoint of oj. 

The dilation ensemble ^ includes SAWs of all lengths. To simulate it using the pivot 
algorithm we must relate it to the ensemble of SAWs with a fixed length in which each SAW 
has the same weight. Our method for doing this is closely related to the method in [1] for 
relating the ensemble of SAWs in the half-plane with a fixed number of steps to radial SLE. 
We decompose the sum in Z(if>), as defined in according to the length of us. 

Z(V0 = J>~ n £ l c (a;)l(A 1 <A(a;)<A 2 )W(a;)V(^y) 

n w.\{jj\=n 
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Let c n be number of SAW in the full plane starting at with n steps. We have 
Z^) = J2cn^ n - J2 ^Ml(Ai<AM<A 2 )%)^(-^-) 

n ui:\w\=n 

The constraint A (a;) > Ai implies that u must have at least Xi/8 steps. So as the lattice spacing 
goes to zero, the first n for which the summand in the sum on n is nonzero goes to infinity. 
Since c n is asymptotic to //"n 7_1 , we replace c n p~ n by n 7-1 . Let P n be the uniform probability 
measure on all SAWs starting at the origin with n steps, and let E n be the associated expected 
value. Then we can write the above as 

Z(V) = Y^tf- 1 E n [\ D (u) l(Ai < A(w) < A 2 ) W(w)V(#t)] 

As noted before, the constraint A (a;) > Ai implies that the sum on n is only over large values. 
So we should be able to approximate E n by its scaling limit. The constraint Id (a;) is a bit 
tricky in the scaling limit, so we proceed as follows. Fix a large positive integer N. If n is also 
large and u and 7 are drawn from P n and P/v, respectively, then the distributions of 8~ 1 n~ u uj 
and 5~ l N~ Ur y are approximately the same. So we will replace E n by En by replacing b~ x n~ v u 
by 8~ 1 N~ v 'y, i.e., we replace co by n v N~ v ^. So A (a;) becomes \(n u N~ Ur y) = n u N~ v X(j) . Our 
weight VK(w) only depends on the polar angle of the endpoint of u>, so we can replace W(cu) by 
Win). And we can replace VCjgy) b Y ^(a^)- 

The constraint Id (a;) is more subtle. The probability that an n-step SAW stays on one side 
of a half plane is conjectured to go to zero as n~ p as n — > 00 with p = 25/64 [7]. So we expect 
that the probability that Id(^) — 1 also goes to zero as n~ p . So we approximate ti p 1d(uj) with 
A^ p 1d(7), i.e., we replace Id{w) by N p n~ p lo{l)- We now have 

Z(V) « iV^n 7 - 1 -^ w [lD(7) l(Ai < n v N- v \{-y) < A 2 ) W( 7 ) V>( w^r)] 

The n dependent part of this is 

^n 7 " 1 -" l(Ax < n v N- v X(i) < A 2 ) 

n 

Since the constraint restricts the sum to large values of n, n' y ~ l ~ p is slowly varying and so we 
can think of this as a Riemann sum approximation to 







x <r-i- P < rjM N -^x{^) < A 2 ) rfx = ciV 7 - p [A( 7 )] (p - 7)/i/ 
where the constant c depends on Ai and A 2 , but nothing else. Thus 

Z(V) ~ ciV 7 ^^[A(7) (p - 7)/ " 1d(7) W( 7 ) ^(ttt)] 

A(7j 
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By taking if) = 1, this result also gives an expression for Z. 

The above derivation was for the radial case in which the starting point of the SAW is in 
the interior of D. If we take a domain which has the origin on its boundary and consider SAW's 
which start at the origin, stay in D and end on its boundary, then we can repeat the above 
derivation. The one important difference is that the number of walks c n should now grow like 
p, n n y ~ 1 ~ p , rather than p^n 1 ^ 1 . So the appropriate power for A (7) becomes p = ^p~- 

We conjecture that when we take the scaling limit our approximations become exact. More 
precisely, we make the following conjecture. 

Conjecture : Let D be a simply connected domain which contains the origin (radial case) 
or has the origin on its boundary (chordal case), and which is star shaped with respect to the 
origin. Let if) be a function on simple curves in D which start at the origin and end on the 
boundary. Let Z(ip)/Z be the expected value of if) in the dilation ensemble of SAW's on a lattice 
of spacing S as defined by Let be the uniform probability measure on N-step SAW's 
7 starting at the origin in the half plane (chordal case) or in the full plane (radial case). Let 
A(7) > be such that 7 ends on X(-j)dD. Let Id (7) be the indicator function of the event that 
7 is inside the domain X(-y)D except for its endpoint(s). Define the weight W{pf) by (T3|). Then 

Zty) E N [\{i?\ D {i)W{i)1>{-fa)\ 
hm — - — - = iim — — ( 4 ) 

i_v z iv->oo E N [x{^yi D {^) w(i)\ 
p — 7 

p = = —61/48 (radial case), 

v 

p = — — - = —3/4 (chordal case). 



3 Lattice effects in boundary densities 

We expect that in the scaling limit the endpoint of the SAW on the boundary of the domain 
D will have a distribution that is absolutely continuous with respect to arc length along the 
boundary. We denote this boundary density by p. If we consider ordinary random walks instead 
of self-avoiding walks, then in the scaling limit this boundary distribution would be harmonic 
measure. 

Let / be a conformal map on D. Then it is expected that the boundary density transforms 

as 

p D (z) = c\f(z)f*p }{D) (f(z)), zedD (5) 

except for a local lattice effect that persists in the scaling limit. If we take g to be the conformal 
map of D onto the unit disc that fixes the origin, then the above implies that Pd(z) is propor- 
tional to \g'(z)\ 5 ^ 8 . Note that the conformal invariance of harmonic measure implies that the 
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boundary density for the ordinary random walk transforms in an analogous way except that 
the power of 5/8 is replaced by 1. Consequently (jSJ) implies that the boundary density for the 
SAW is proportional to the boundary density for the ordinary random walk raised to the 5/8 
power. 

We now turn to the computation of the lattice effect correction to the boundary density. 
This computation for the dilation ensemble is analogous to what was done in [5]. The constraint 
that the SAW u stays inside the dilated curve X(w)dD has both a macroscopic and microscopic 
nature. The conjecture (jSJ) is the result of the macroscopic effect. Near the endpoint of the 
walk there is an additional microscopic effect. Consider a SAW that ends at z and consider the 
tangent line to the curve 3D at z. The constraint that the SAW stays inside X(u)D implies 
that near the endpoint the SAW must stay on one side of this line. This will produce a factor 
1(9) that depends on the angle of the tangent line with respect to the lattice. 

For an angle 9, let L be the line with polar angle 9 passing through the origin. We consider 
walks with N steps starting at the origin. Let be the number of such walks, and let b^(9) 
be the number of such walks that stay on one side of the line. So b^(9)/cN is the probability 
that an N step SAW stays on one side of the line. We expect that this probability goes to zero 
as N~ p as N — > oo, and we conjecture that the lattice effect function is given by 

m = Km h J^ N P ( 6 ) 

Since 9 and 9 + 180 give the same line, 1(9) has period 180 degrees. The function will have 
more symmetries depending on the type of lattice. For example, for the square lattice 1(9) has 
period 90 degrees and 1(9) = 1(90 — 9). If we use the weight W(u) in ([1]) and (T5]), then the 
boundary density for the dilation ensemble will be pn(z)l(9(u)) where 9(u) the polar angle 
of the endpoint of u>. To remove this lattice effect from our dilation ensemble we replace the 
weight W(uj) given by ([3]) by 

In our conjecture (jlj) the weight W(u) is used in both sides of this equation. The same derivation 
shows that this equality should also hold if we use W(u) in both sides. 

We take the constraint that u> stays inside X(u)D to mean that it stays strictly inside the 
curve except for the endpoint. The other convention would be to allow u> to have sites that 
lie on X(u)dD in addition to its endpoint. If 3D has flat sections, then there can be a big 
difference between these two conventions for angles 9 such that a line through the origin at 
angle 9 passes through some lattice sites. This will be the case for the equilateral triangle that 
we consider in our simulations. 

In figured] we plot the lattice effect function 1(9) for the dilation ensemble and for two other 
ensembles, all on the square lattice. One is the cut-curve ensemble studied in [5]. The other 
ensemble is the "natural" ensemble in which we take all SAWs such that all the sites on the 
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Figure 1: The lattice effect function 1(9) for three different ensembles on the square lattice. 



walk are in D except for the endpoint which is outside of D. In other words, the last bond of the 
walk crosses the boundary of D and is the only bond which does so. All three functions have 
been normalized so that the total area under each curve is 1. It is important to note that the 
vertical scale of the plot does not include 0. These functions are actually relatively flat. Note 
that the lattice effect function for the dilation ensemble studied in this paper is considerably 
flatter than the function for the cut-curve ensemble that was simulated in [5]. 

The function 1(9) is continuous for the cut-curve and natural ensembles, but for the dilation 
ensemble it is not. For the square lattice it is discontinuous at 9 such that tan(#) is rational. 
The biggest discontinuity is at 9 = 0, but it cannot be seen in the figure since 1(0) = 0.001516 
which is well below the region plotted. The second biggest discontinuity at 9 = 45 is seen 
in the figure as an isolated point at /(45) = 0.002361. This discontinuity is atypical in that 
1(4:5—) = Z(45+). Other discontinuities with 1(9— ) ^ 1(9+) can be seen in the figure. To see 
why 1(9) is discontinuous, consider the case of 9 = 0. Recall that we make the convention that 
the walk must stay strictly inside the domain except for the endpoint. So when 9 is exactly 0, 
b^(9) counts walks that start at the origin and then stay strictly on one side of the horizontal 
axis. So the walk cannot visit any site on the horizontal axis. Now consider a 9 that is slightly 
greater than zero and consider a walk which stays above this line. The sites on the negative 
horizontal axis are now above the line and so the walk can visit these sites. Thus b^(0) is 
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significantly smaller than &tv(0 + ). If we change the convention that the SAW must stay strictly 
inside the domain to allow SAW's that visit sites on the boundary, then the values at the angles 
where 1(9) is discontinuous will change but it will still be discontinuous. 

4 Simulations 

We use our conjecture to simulate the dilation ensemble and compare the boundary density 
found in the simulation with the density given by (j5J). Note that our simulations are testing three 
different conjectures. One is the conjecture (Hj) that says we can use the fixed length ensemble 
to simulate the ensemble of SAWs in a bounded domain. Another is the SLE partition function 
prediction for the boundary density (j3j). And finally there is the conjecture that the lattice 
effect correction to this density is given by 1(9), i.e., by using the weight (|7j). 

The exact predictions for the boundary density from SLE partition functions (jSJ) are given in 
the appendix. In our simulations we work with cumulative distribution functions (cdf 's) instead 
of densities. Extracting a density from the simulation requires taking a numerical derivative, 
i.e., choosing a bin size and computing a histogram. Using the cdf's avoids this. 

In our conjecture (jl]) the probability measure En only depends on whether the geometry 
we are studying is chordal or radial. So we only need to do two simulations, one for the chordal 
cases and one for the radial cases. In each simulation we sample the Markov chain every 100 
time steps. In the chordal case we generated 9 x 10 9 samples and in the radial case 11 x 10 9 
samples. For each sample and each domain D we test if 1d(w) = 1. When it does we use that 
sample in the computation of the cdf of the boundary density for that domain. The probability 
that 1d(oj) = 1 depends on N and on the domain. In our simulations it ranges from |% to 2%. 

The first domain we consider is a horizontal strip of height 1. For this domain there are no 
lattice effects since the entire boundary is horizontal. We consider both a chordal case and a 
radial case. In the chordal case the strip is {z : < Im(z) < 1}. The SAW starts at the origin 
and ends on the upper boundary. In the radial case the strip is {z : —1/4 < Im(z) < 3/4}. 
The SAW starts at the origin and can end on either boundary. The conjectured density for 
the chordal case was tested by simulation in [3|. Here we are primarily interested in using 
these domains to test our conjecture (jl]) that uses the fixed length ensemble for the simulation, 
especially the value of the power p. 

In figure [2] we plot six curves. Two curves are the cdf's computed using SLE partition 
functions for the chordal and radial cases. The other four curves are simulation cdf's for the 
chordal and radial cases. In each case we have two simulations, one with the power p given by 
flSJ) and one with p — 0. The curves from the SLE prediction and the curves from simulations 
with the correct value of p are indistinguishable on this plot, and so it appears there are only 
four curves in the figure. The maximum difference between these curves is given in table [H It 
is on the order of 10 -4 in the chordal case and 1.5 x 10~ 4 in the radial case. The discrepancy 
for the curves with p = is quite large for the radial strip. For the chordal strip it is smaller, 
but can still be clearly seen in the plot. The maxima of the differences for the p = case is 
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Domain 




1(9) used ? 


max of difference of cdf's 
(in thousandths) 


St.Tin wit.n t) = 


rhnrrlal 




26 90980 ± 63765 


Strip with p = 


radial 




229.35605 ± 0.64829 


Sif TlFl 

k-> LI 1 kJ 


priArn fil 




09842 + 1 1 61 9 


Strip 


radial 




15563 ± 15085 

V 7 • _L k./ k./ V 7 v7 — 1 — V 7 • -L- v7 V 7 V 7 '.7 


J. 1 idiitiic 


T*S1 H 1 £1 1 
1 CXU.1C11 


yes 


05939 4- 08689 


Centered circle 


radial 


no 


3.16364 ± 0.07145 


Centered circle 


radial 


yes 


0.10997 ±0.07793 


Off-center circle 


radial 


no 


2.26140 ±0.09184 


Off-center circle 


radial 


yes 


0.07513 ±0.08061 


Partial circle 


chordal 


no 


1.28211 ±0.06962 


Partial circle 


chordal 


yes 


0.09512 ±0.06158 


Tangent circle 


chordal 


no 


1.87370 ±0.10155 


Tangent circle 


chordal 


yes 


0.18705 ±0.10116 



Table 1 : For each domain the last column is the maximum of the absolute value of the 
difference between two cdf's. One cdf is computed using the SLE partition function 
prediction for the density. The other is from the SAW simulation. The middle column 
indicates if the lattice effect 1(9) is corrected for. 



also given in the table. 

We can estimate the power p by minimizing the difference between the boundary density 
computed by simulation using the conjecture (jlj) and the SLE prediction for the boundary 
density. We use the two strip geometries described above and minimize the L 2 norm of the 
difference between the densities. Our simulations to estimate the power p are separate from 
the simulations to compute the boundary densities. We use SAWs with 100, 000 steps. For the 
chordal geometry we generated approximately 8 million samples and for the radial geometry 
approximately 10 million samples. By samples we mean SAWs for which 1d(^->) = 1- lu the 
chordal case the conjectured exact value is p = —3/4 = —0.75 and from the simulations we 
find the minimum difference is when p = —0.751874, a difference of 0.25%. In the radial case 
the conjectured exact value is p = —61/48 = —1.270833 and from the simulations we find 
-1.269917, a difference of 0.07%. 

Next we consider an equilateral triangle centered at the origin whose vertices have polar 
angles of 0, 120 and 240. The side corresponding to polar angles in [120, 240] is vertical and 
the sides corresponding to [0, 120] and [240, 360] are at 30 degrees with respect to the lattice 
directions. So two of the sides will have the same value of 1(9) while the third side has a 
different value. This geometry gives an extreme example of the lattice effects. We simulate 
the ensemble with the weight factor W(u) which does not correct for the lattice effects. We 
find that the probabilities for hitting the sides corresponding to [0, 120], [120, 240], [240, 360] are 
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Figure 2: Comparison of cdf 's from simulations and theory for chordal strip (polar 
angle ranges over [0, 180]) and radial strip (polar angle ranges over [0, 360]). The solid 
curves are simulations using p = 0. The overlying dashed curves are the simulations 
with the correct p and the SLE prediction. 

0.387375, 0.225173, 0.387452 respectively. The ratio of the smaller probability to the average 
of the other two probabilities is 0.581221. This should be compared with the ratio of the two 
values of the lattice effect function which is Z(0)//(30) = 0.581281. In figure [3] we plot the cdf 
from our simulation. The smaller probability of hitting the vertical side is clearly seen. The two 
horizontal lines are at the heights predicted by the lattice effect function. For this geometry 
the only lattice effect should be to make the probabilities of hitting the three different sides 
unequal. Given that you hit a particular edge, the distribution along that edge should be the 
same for the three edges. So we can remove the lattice effect by looking at the polar angle of 
the endpoint mod 120. We compute the cdf of this random variable. In the inset in figure |3] we 
show this cdf minus the SLE partition function prediction for this cdf. This difference is quite 
small, on the order of 5 x 10~ 5 . The maximum of the difference is given in tabled] 

The next domain we consider is a unit disc centered at the origin. As always, the SAW starts 
at the origin. The SLE partition function prediction is just that the density is uniform. The 
deviation from this is entirely due to lattice effects. The boundary density will have a period of 
90 degrees, so we take the polar angle of the endpoint of the walk mod 90. In figure @] we plot 
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Figure 3: Simulation cdf for the equilateral triangle. The inset shows the difference 
of the simulation cdf and the theoretical cdf when the polar angle is reduced mod 120 
to remove the lattice effects. 

two curves. Both are the cdf from the simulation minus the cdf of the uniform distribution. 
In one we use the weight factor W(co) which does not include the lattice effect correction, and 
in the other we use the weight factor W(u>) which does include the correction. The sine-like 
curve is for the simulation without the lattice correction. The difference is small, on the order 
of 0.003, but clearly not zero. The flat curve is the difference when we do include the lattice 
effect correction in the simulation. The size of this difference is a test of our prediction for the 
lattice effect correction. The difference is extremely small, on the order of 10~ 4 . The maximum 
of the difference with and without the lattice correction is given in table [U 

Our last three simulations involve geometries for which the conjectured density from (jSJ) 
is non-trivial and the lattice function 1(9) enters in a non-trivial way. In the first of these 
geometries the domain is a unit disc centered at 3/4. We consider SAW's that start at the 
origin and end on the arc of the circle. As before we consider the difference between the cdf 
from the simulation and the cdf predicted by SLE partition functions. As before we consider 
two cases, one using the weight W(oS) and one using the weight W(u). We refer to this geometry 
as the "off-center circle" and give the maximum of the difference of the cdf 's in table [TJ As 
in the previous geometry the difference when we do not correct for the lattice effect is clearly 
non-zero, while it is zero within the errors in the simulation when we do correct for the lattice 
effects. 
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Figure 4: Difference of cdf's from simulations and theory for disc with origin at the 
center. The sine-like curve is for the simulation which does not correct for the lattice 
effects; the flat curve is for the simulation that does correct. 

In the next geometry the domain is the intersection of the unit disc centered at — 3i/4 and 
the upper half plane. So the domain is bounded by the real axis and the portion of the unit 
circle centered at — 3z/4 that lies above the real axis. We consider SAW's that start at the 
origin and end on the arc of the circle. This geometry is called "partial circle" in the table. 
In the final geometry the domain is a unit disc centered at i. So the real axis is tangent to 
the disc at the origin. We consider SAWs starting at the origin and ending on the upper half 
of the boundary. So the polar angle with respect to the origin ranges from 45 to 135 degrees. 
This geometry is called "tangent circle" in the table. For both of these geometries the maxima 
of the differences is given in table [TJ They are clearly non-zero when we do not correct for the 
lattice effects, while they are zero within the errors in the simulation when we do correct for 
the lattice effects. 

5 Conclusions 

We have studied SAWs in a bounded domain D which start at a fixed point (either in the 
interior or on the boundary) and end at an unconstrained point on the boundary. The dilation 
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ensemble of SAW's consists of all SAW's starting at the origin with the property that they can 
be scaled to give a curve inside our domain D which ends on the boundary of D. It can be 
thought of as interpreting the constraint that a SAW ends on the boundary of D by thickening 
the boundary. We weight a SAW u by the usual weight of times two factors. The factor 
W(u), given by (El), accounts for the varying thickness of the boundary. The other factor of 
1/1(0) accounts for the local lattice effect near the endpoint of the SAW that persists in the 
scaling limit. 

We have conjectured that the scaling limit of the dilation ensemble is related to the limit as 
iV — > oo of the uniform measure on SAW's with iV-steps conditioned on the event that when the 
SAW is scaled so that it ends on the boundary of D, the SAW lies entirely in D. We have used 
our conjecture to simulate the SAW in several bounded domains and compared the distribution 
of the endpoint of the SAW on the boundary with the SLE partition function prediction of 
this density. We find excellent agreement. This supports our conjectured relationship between 
the dilation ensemble and the uniform probability measure on SAW's of a fixed length, the 
conjecture for the lattice effects that persist in the scaling limit and the SLE partition function 
predictions for the boundary density. Note our simulations only looked at the boundary density 
of the endpoint of the SAW. By using the SLE prediction for the boundary density and various 
explicit predictions about chordal and radial SLE, one can make explicit predictions about 
the SAW curve inside D and use our conjecture to test them. We have not carried out such 
simulations. 



A Computation of SLE partition functions 

In this appendix we give the analytic computations of the boundary density for the various 
domains. We use p(-) to denote the boundary density with respect to arc length (or a constant 
multiple of arc length) along C . In all cases the SAW starts at the origin. The constant c in 
the following is determined by the requirement that p is a probability density. 

A.l Strip - chordal 

The domain is a strip of height 1 with the origin on the lower boundary. We are interested in 
the ensemble of SAW's in the strip that start at the origin and end at some point on the upper 
boundary of the strip. We parametrize the upper boundary by x + i where — oo < x < oo. A 
straightforward application of (jSJ) yields the boundary density with respect to dx. 

p(x) = c [cosh(7ra;/2)]~ 5 / 4 

A. 2 Strip - radial 

The domain is a horizontal strip of height 1, with the origin at a height h above the bottom 
boundary. So the domain is {z : —h < Im(z) < 1 — h}. We consider the ensemble of SAW's in 
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this strip that start at the origin, and end at a point on either the top or bottom boundary of 
the strip. We parametrize the lower boundary by x — ih and the upper boundary by x + (1 — h)i 
where — oo < x < oo. Then the boundary density with respect to dx is 



c [cosh(7rx) — cos(nh)] 

C [cOsh(7Tx) + COs(7l7l)] 



5/8 



if y = -h 

if y = l-h 



p(x + iy) 



5/8 



A. 3 Triangle 



We use the Schwarz-Christoffel mapping 



F(z) = (w + 1)" 2/3 (w - l)- 2/3 dw 



■1 z 



which maps the upper half plane onto an equilateral triangle. It sends —1, 1, oo to the vertices, 
and some simple substitutions show it sends —3, 0, 3 to the midpoints of the three sides. To 
find the pre-image of the center of the triangle, note that the center of triangle is the unique 
point fixed by rotations about the center. These rotations are the conformal automorphisms of 
the triangle that permute the vertices. So they correspond to conformal automorphisms of the 
half plane that permute —1, 1 and oo. Such Mobius transformations fix iv3, so F(iy3) is the 
center of the triangle. 

A simple application of (jSJ) shows that for SAWs in the half plane starting at iy?>, the 
boundary density along the real axis is proportional to (x 2 + 3) _5//8 . By symmetry it is enough 
to find the boundary density for the triangle between the midpoint of one edge and a vertex. 
We consider F(3, oo). (For \z\ > 3, F(z) can be computed numerically very quickly by a power 
series expansion.) Let Iq be the length of the edges in the triangle. We define I — l(x) by 



Then l(x) is the distance from F(x) to the vertex F(oo). Applying fl5]) we find the boundary 
density with respect to dl: 



A. 4 Off-center circle 

The domain D is a unit disc centered at a + ib where |a + zfe| < 1 so that the origin is inside the 
disc. Let <fi be the polar angle of a point on the boundary with respect to the center at a + ib. 
So 4> is proportional to arc length along the boundary. An application of (JSJ) with a Moibius 
transformation yields the boundary density with respect to d<p. 



2Z(» F(x) 



l ~ F(3) 



p(l) = c 



| F /( x )|-5/8 ^ - 1)5/12 

(x 2 + 3) 5 / 8 ~ C (x 2 + 3) 5 / 8 



p(4>) = c [l + a 2 + b 2 + 2a cos (j) + 2b sin 



5/8 
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A. 5 Partial circle - chordal 



We consider the unit disc centered at b where b is real and |6| < 1. So part of the disc lies 
below the real axis. We take the domain to be the intersection of this disc with the upper half 
plane. We consider SAWs in this domain that start at the origin and end on the arc of the 
circle above the real axis. 

The circle intersects the real axis at ±d with d = yl — b 2 . The map — (z + d)/(z — d) maps 
the domain to the wedge < arg(z) < (3 where tan/3 = — yl — b 2 /b. So if we let 



/(*) 



z + d 
z — d 



7T//3 



then / maps the domain to the upper half plane and sends to 1. The endpoint of the walk is 
mapped to a point on the negative real axis. 

An easy application of (JSJ) shows that for the SAW in the upper half plane starting at 1 and 
ending at x on the negative real axis, the boundary density is proportional to (1 — a;) -5 / 4 . So 
applying (JSJ) to the map / shows that the boundary density with respect to d<p is 



p{4>) = c 



\z + d\ 



5/8 



Jz-dlWJ (l-f(z)) 5 / 4 



where z = ib + e 1 ^. 



A. 6 Tangent circle - chordal 

The domain is a disc of radius 1 centered at % so it is tangent to the real axis at the origin. Let 
<j) be the polar angle with respect to i, so is proportional to arc length along the boundary. 
We condition on the event that the walk stays in this disc and ends on the upper half of the 
circle bounding the disc, i.e., the arc of the boundary where < <fi < 180. 

The boundary density for the upper half plane when the SAW starts at the origin is x _5//4 . 
So using a Mobius transformation to map the disc to the half plane, the formula (jSJ) gives the 
density with respect to d(f>: 

p(0) = c \l - cos0 + sin0|- 5/4 (l - cos0) 5/8 
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